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Abstract. We present a CSO and BIMA molecular line 
survey of the dense, quiescent molecular environment 
ahead of HH 2. The molecular gas is cold, 13 K, and 
moderately dense, 3 10 5 cm -3 . A total of 14 species has 
been detected (including different isotopes and deuterated 
species). The relative abundances of the clump are com- 
pared with other dense molecular environments, includ- 
ing quiescent dark clouds, and active low and high mass 
star forming regions. This comparison confirms the pe- 
culiar chemical composition of the quiescent gas irradi- 
ated by the HH objects. Thus, from this comparison, we 
found that the HCO+, CH 3 OH and H 2 CO are strongly en- 
hanced. SO and SO2 are weakly enhanced, whereas HCN 
and CS are under abundant. The CN abundance is within 
the range of value found in starless dark clouds, but it 
is low with respect to high mass star forming regions. 
Finally, the chemical composition of HH 2 confirms the 
qualitative results of the Viti & Williams (1999) complex 
chemical model that follows the chemical behavior of a 
molecular clump irradiated by a HH object. 
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1. Introduction 

Herbig-Haro (HH) objects, which trace shock-excited 
plasma, are signposts of the intense outflow phenomena 
associated with star formation (e.g., Reipurth & Bally 
2001). Several HH objects are found to have associated 
quiescent dense clumps ahead of them: HH 1/2 (see ref- 
erence s belo w), HH 7-11 (R udolph & Welch |1988[ Dent 
et al. |j~993;; Rudolph et al. pOOl] ), HH 34 (R udolph & 
Welch p92| , NGC 2264G (Girart et al. [2000| ), HH80N 
(Girart et al. |1994[ |200lt Girart, Estalella & Ho |199Sp . 
Most of them have high excitation knots, i. e., with strong 
high excitation lines, such as [OIII] and Ha, strong 
UV lines and significant UV continuum radiation (e.g., 



Reipurth & Bally 2001). The quiescent clumps are charac- 



terized by the low temperature and narrow line widths of 
the molecular emission, and by the emission enhancement 
of HCO + and NH3. In spite of the presence of the HH ob- 
jects, there is a lack of apparent dynamical perturbation in 
the clumps. However, recent observations show clear sign- 
posts of sta r formation within the HH 80N clumps (Girart 
et al. |2001| ). 



Girart et al. (1994) suggested that the regions of en- 
hanced HCO + and NH3 emission were a consequence of 
the irradiation from the HH shock affecting the chemistry 
within a small dense clump nearby in the molecular cloud. 
It was supposed that the UV radiation would evaporate 
icy mantles and promote photochemistry in the enriched 
gas. Taylor & Williams ( 199(f ) modeled this situation with 
a simple chemistry and showed that HCO + would arise 
from the interaction of C + (derived from CO) with H2O 
liberated from the ice. The chemical enhancements would 
be transient, but predicted column densities were signifi- 
cant. Viti & Williams ( |1999| , hereafter VW99) noted that 
the impact of a high radiation field on a dense gas should 
lead to a rich chemistry, and they therefore extended the 
1996 model. They predicted that, in addition to HCO + 
and NH3, a wide variety of species should arise with en- 
hanced abundances, including CH3OH, H2CO, SO, SO2, 
and CN and that the effects in a single clump should last 



for 10 4 years. Raga & Williams (2000) showed that the 



photochemical effects created by the HH object should re- 
spond to its movement through the molecular cloud, and 
predicted that the enhanced regions should have a char- 
acteristic morphology. 

Implicit in the model is the idea that molecular clouds 
are clumpy on a small scale. This was suggested by Tay- 
lor, Morata & Williams ( |1996[ ) as a necessary assump- 
tion to interpret the relative spatial distributions of CS 
and NH3 in molecular clouds in terms of time-dependent 
chemistry. The clumps were required to be transient on 
a timescale of about one million years. Subsequent inter- 
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ferometer observations (e.g. Peng et al. 1998 , for TMC-1 
Core D; Morata, Girart & Estalella |2002| for L673) have 
supported the Taylor et al. model. Therefore, the HH ob- 
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Table 1. Lines Observed with the 10-m CSO telescope 
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a Av is the FWHM line width 

b Overlapped with the NS 11/2,11/2-9/2,9/2 transition 



jects may be regarded as probes of the transient substruc- 
ture of molecular clouds. Wherever the probe encounters 
a small clump, a characteristic chemistry arises. This pa- 
per reports a search for that characteristic chemistry in 
the clump associated with HH 2. 



The prototypical HH 1 and HH 2 objects (Herbig 1951 



Haro 1952), located in the L1641N molecular cloud in 
Orion, provided the earliest spectral evidence that the HH 
objects trace strong shocks due to the presence of stellar 
winds (Schwartz |1978| ). Both HH 1 and HH 2 belong to 
the class of high excitation HH objects, with strong UV 
(Bohm ct al. 1987; Raymond, Hartmann & Hartigan 1988; 
Raymond, Blair & Long 1997) and centimeter continuum 
emission (Rodriguez et al. |1990| , |2000| ). HH 2 shows a very 
complex and apparently chaotic morphological, kinemat- 



ical and excita tion s tructure (e.g., Schwartz et al. 1993; 
Eisloffcl et al. 1994 ). Detailed HST observations reveal 



that this complex structure is the consequence of the in- 
teraction of the HH 2 jet with a dense ambient cloud, 
resulting in a very bright, high ionization and complex 
structure (Hester, Stapelfeldt, & Scowen 199S). Recent 



Chandra observations reveal that the strongest HH 2 knot, 
which is also the highest excitation knot, has associated 
X-ray emission (Pravdo et al. 2001). 



The presence of quiescent dense ambient gas ahead 
of HH 2 has been well established from HCO + and NH 3 



observations (Davis, Dent & Bell Burnell 199C ; Torrelles et 
al. |1992| , |1994| ; Choi & Zhou |l997[ ). The ammonia emission 
shows a clumpy medium, with clump sizes of only 20" 



(Torrelles ct al. 1992) or ~ 8000 AU (assuming a distance 
of 390 pc: Anthony- Twarog [19821). JCMT observations by 



Dent (1997) shows that the chemical composition of the 
dense ambient gas is strongly altered by the presence of the 
HH 2 radiation, with a str ong em ission enhancement of the 
HCO + . Wolhre & Konigl ( 1993 ) carried out a shock model 
where the X-ray, EUV and FUV radiation field from HH 2 
penetrates the cloud, initiates an ion chemistry, enhancing 
the abundances of the HCO + and the electrons, which 
produces an excitation enhancement of the HCO + . The 
combination of these two effects lead to a strong emission 
enhancement. 

In this paper we present 0.8-1.4 mm CSO observa- 
tions towards the quiescent clump ahead of HH 2. We also 
present 3 mm spectra obtained with the BIMA array. The 
goal of this paper is to study the chemical effect produced 
by the strong HH 2 radiation over the dense quiescent am- 
bient gas around HH 2. A second paper will analyze and 
discuss the spatial distribution of the molecules observed 



with the BIMA array (Girart et al. 2002, hereafter Paper 
II). 
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Table 2. Lines Observed with the BIMA Array 
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a Av is the FWHM line width 



2. Observations and Results 

Observations were carried out in October 1999 with the 
10.4 m telescope of the Caltech Submillimcter Observa- 
toryQ (CSO). The 200-300 GHz and 300-400 GHz re- 
ceivers were used in conjunction with 50 MHz and 500 
MHz acousto-optical spectrometers, which provided ve- 



1 The Caltech Submillimeter Observatory is funded by the 
National Science Foundation under contract AST-9615025 



locity coverage of 60 and 600 km s^ 1 at 250 GHz and 
of 40 and 400 km s _1 at 350 GHz. Pointing, focus and 
calibration were checked every few hours using Saturn. 
The main-beam efficiency, measured by observing Sat- 
urn, were -67% at 250 GHz and -60% at 350 GHz and 
system temperatures ranged 240-670 K. All the spectra 
were observed at the position: a(J2000) = 5 h 36 m 26^84; 
<5(J2000) = -6°47'27'.'l. This position is very close (~ 3") 
to the HH 2 knot L (e.g., Hester et al. |99f) and to the 
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Fig. 1. Spectra of the detected lines with the BIMA and CSO telescopes. Lines from BIMA were convolved with a 
Gaussian in order to obtain a resulting beam of 30", approximately the beam of the CSO data at 250 GHz. 



peak emission of the HCO + 4 3 integrated map from Dent 
( |1997| ) . Table [l] lists all the molecular lines observed along 
with their rest frequency, the angular resolution at the 
observed frequency and the Gaussian fits to the line pro- 
files. The upper limits show the 3-<r level within the 5.2- 
7.8 km s _1 range. 

The 10-antenna BIMA arrayQ observations were car- 
ried out between 1999 October and 2001 May in the C 
configuration. A detailed description of the observations 



2 The BIMA array is operated by the Berkeley-Illinois- 
Maryland Association with support from the National Science 
Foundation. 



and reduction is given in Paper II. In order to compare 
with the CSO spectra, we smoothed the BIMA maps by 
convolving them with a Gaussian, with a resulting FWHM 
beam of 30", the angular resolution of the CSO observa- 
tions at 230 GHz. A spectrum of each transition observed 
with BIMA was synthesized at the position observed with 
the CSO telescope (see paragraph above). Table || lists all 
the molecular lines observed, their rest frequency, and the 
Gaussian fits to the line profiles. 

Figure |l] shows the spectra of the detected lines with 
the CSO and the BIMA observations smoothed to the 
30" angular resolution. A total of 33 transitions (including 
hyperfine transitions) were detected for 14 species: 13 CO 
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5 



and C 18 0, HCO+, H 13 CO+, HC 18 0+, HCN, H 2 CO, CN, 
CH3OH, CS, SO, S0 2 and the deuterated species DCO+ 
and DCN. The spectral lines peak at a vlsr within the 
range of - 6.4 km s" 1 (CS and HCN) to ~ 6.9 km s" 1 
(HCO + ). This range of values is in agreement with pre- 
vious observations of the quiescent gas ahead of HH 2 



-6°46' 



(Torrelles et al. [199$ Choi & Zhou |l997| ). The line widths 
range from ~ 0.8 km s _1 for the DCO + to ~ 1.9 km s _1 
for HCO+ and HCN. The 1 mm CS and HCN transitions 
have large line widths, ^2.4 km s _1 . However, this value 
is not significant due to the low SNR of these lines. The 
differences in line center velocity and line width are signif- 
icant and are due to the combination of velocity gradients 
within the quiescent ambient gas and a spatial chemical 
differentiation, which will be discussed in Paper II. Yet, 
for the HCO + the larger line width can also be due to a 
high optical depth and the presence of high velocity emis- 
sion. The HCO + 1-0 spectrum shows clearly rcdshifted 
and blueshifted wings, which implies the presence of high 
velocity molecular gas. This high velocity component is 
tracing the interaction between the HH object and the 
clump, which confirms the spatial association of the clump 
with HH 2 (Paper II). 

2.1. HH2 and the clump 

Figure |^ shows the distribution of the quiescent molecular 
gas, traced by the 13 CO emission, ahead of the HH 2 re- 
gion. There is clear column density enhancement, mainly 
downstream of HH2, whereas at the position of HH 2 or 
just behind it there is little ambient gas. The SO emission 
is also shown to better locate the chemically enhanced 
region just ahead of HH 2, alm ost coincident with HH 2 
knot L (e. g., Hester et al. |l998[ ). The high excitation HH 2 
emission, as well as the X-ray and the radio continuum 
emission, arise from the brightest and largest knot in the 
[SII] image (see Fig|). 

3. Analysis 

3.1. The Physical Conditions of the Quiescent Gas 

HCO + , SO and CH3OH are the molecules with more than 
three transitions detected (including the rarer isotopes for 
the HCO + ). Therefore, the observed emission from these 
molecules can be used to derive the physical condition 
of the quiescent gas. SO and CH3OH are molecules with 
moderate dipole moments, 1.55 and 0.90 Debyes, respec- 
tively, which implies that if the densities are not too low 
their emission will be close to thermalization. Thus, we 
carried out the population diagram analysis (e.g., Gold- 
smith & Langer 1999) to derive the temperature of the 



gas and the column densities of SO and CH3OH. For the 
HCO + , since we have two rare isotopes observed, H 13 CO + 
and HC 18 + , we carried out the Monte Carlo radiative 
t ransf er model developed by Hogerheijde & van der Tak 
(|2000D . 



-6 U 47' 



-6"48' 




5 h 36 m 30?0 



25 s 



Fig. 2. Superposition of the gray scale of the [SII] im- 
age (from Curiel, private communication), the integrated 
emission of the 13 CO 1-0 line (thin contours) over the 5.1- 
8.1 km s _1 «lsr interval, and the integrated emission of 
the SO 32~2i line (thick contour) over the aforementioned 
«lsr range. Contours levels are 4, 13, 25, 37, 49, ... times 
the rms noise of the map, 47 mJy beam^ 1 for the 13 CO, 
and 50% of the peak intensity for the SO. The synthesized 
beam of the 13 CO (left) and the SO (right) are shown in 
the lower left corner. The star marks the position of HH 2 
VLA 1, the powering source of the HH 1/2 outflow. The 
cross marks the center coordinates of the single-dish spec- 
tra. The two concentric dashed circles show the 21" and 
28" beam sizes of the CSO spectra at 265 and 360 GHz, 
respectively. 

3.1.1. Population Diagram Analysis 

The method used to carry out the population diagram 
analysis, which takes into account approximately the opti- 
cal depth, is described in appendix § |A|. In order to obtain 
more accurate fits for the SO and CH3OH, which were 
observed with two different telescopes, a 20% calibration 
uncertainty was included in the estimated errors. Figure |^ 
shows the population diagram for the SO and CH3OH 
molecules. 

The SO emission is best fitted for T rot ~ 9.4 ± 1.3 K, 



7V(SO)c 



o 9+1.4 
°- z -1.2 



x 10 cm and a filling factor of 



/ = 0.35. This solution gives moderate optical depths 
for the observed transitions, with a maximum value of 
1.0 for the 32~2i transition. The high angular resolution 
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SO 



40 
E, (K) 



Fig. 3. Population diagrams for the SO (left panel) and 
the CH3OH (right panel). Solid filled circles indicate the 
observed data. The error bars include the measurement 
uncertainty and the calibration uncertainty (~ 20%). The 
solid lines show the best fit obtained by assuming opti- 
cally thin emission (i.e., using the standard population 
diagram analysis). The open circles shows the best solu- 
tion for the least square fit taking into account the optical 
depth correction (see § |A|). 




_l_ 
20 



1 

40 

E (K) 



Observed 
Model 



5 10 
Velocity (km/s) 



5 10 
Velocity (km/s) 



Fig. 4. Observed (solid contours) and synthetic (dashed 
contours) HCO+, H 13 CO+ and HC 18 0+ spectra. 



(#fwhm ~ 9") BIMA SO 3 2 -2i maps show that the in- 
tegrated emission arises from a region with a deconvolved 
size of 26" x 15", PA = -34° (Paper II), which implies 
a filling factor of approximately /so — (26 x 15)/(26 x 
15 + 30 x 30) = 0.30. Therefore, the best solution for the 
SO obtained using the method described in appendix § [a], 
gives a filling factor in agreement with the observed value 
from the BIMA maps. 

The population diagram technique for the CH3OH 
emission gives the best fit for T ro t — 10.2 ± 0.9 K, 



A(CH 3 OH)c 



6A±\% x 10 13 



cm" 2 , and / = 0.20. The 



optical depths of the observed transitions are small: the 
maximum value is 0.56, for the 2o-lo A + transition. The 
smaller filling factor of the CH3OH solution with respect 
to the SO is because of its slightly smaller emitting area: 
/ch 3 oh — 0.24 (Paper II). We note that although the op- 
tical depths of the SO and CH3OH are small (5s 1) the 
population diagram technique used here and described in 
§ gives column densities of 20% (CH 3 OH) to 30% (SO) 
higher than the standard population diagram technique 
where optically thin emission is assumed. 

3.1.2. Monte Carlo Radiative Transfer Code 

We used the one dimensional version of the Monte Carlo 
code developed by Hogerheijde & van der Tak ( p000| ) , 
which calculates the radiative transfer and excitation of 
molecular lines. The code is formulated from the viewpoint 
of cells rather than photons, which allows the separation 
of local and external contributions of the radiation field. 



This gives an accurate and fast performance even for high 
opacities (for more details see Hogerheijde & van der Tak 

pooo| ). 

The HCO + emission was assumed to arise from a 
sphere of radius 1 10 17 cm or 6680 AU (17" at the dis- 
tance of HH 2), which is approximately the radius of the 
BIMA HCO + 1-0 emission for the peak intensity compo- 
nent (the HCO + emission extends over a larger region to 
the NE of the peak intensity, but this extended component 
is not picked up by the CSO beam). The volume density 
and temperature were assumed constant and no velocity 
gradient was adopted (we note that the results are insen- 
sitive for radial velocity gradients of ~ 0.3 km s _1 within 
the adopted radius, i. e., ~9 km s -1 pc -1 ). The 12 C 16 to 
13 C 16 and 12 C 16 to 12 C 18 abundance ratios adopted 
were 63 and 560, derived for the Orion region by Langer 



Penzias (1990, 1993] ). We explored a large range of val- 
ues for the volume density, the temperature, the molecu- 
lar fractional abundance and the intrinsic line width. The 
resulting synthetic maps of the different transitions were 
convolved with a Gaussian to match the angular resolu- 
tion of the spectra (28" for the 1-0 and 3-2 transitions 
and 21" for the 4-3 transition). The best fit is obtained 
for a volume density of ~ 2.8 10 5 cm -3 , a temperature of 
~ 13 K, an intrinsic line width of ~ 0.6 km s _1 , and a 
HCO + fractional abundance of 2.0 10 -9 . Figure ^ shows 
the overlap of the observed and synthetic spectra for the 
HCO+, H 13 CO+ and HC 18 0+ transitions. The synthetic 
spectra fit well the observed spectra, except for the HCO + 
1-0. The lack of observed emission at the central and red- 
shifted sections of the line, as compared with the predicted 
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spectra, is possibly due to a combination of two effects. 
First, it could be due to the missing flux, a consequence 
of the lack of short spacing visibilities at the interferom- 
eter. Second, it is possible that the HCO + 1-0 is affected 
by the absorption of a cold and low density component 
of the cloud, whose effects will only be observable in the 
lowest rotational transition of abundant species with high 



dipole moment (see Girart et al. 200C) 



The assumed radius, derived from the HCO + 1-0 
maps, could be highly uncertain due to high optical depth, 
absorption by the foreground cold gas or resolving out the 
large-scale HCO + emission. However, the models run with 
radii higher or lower than roughly a 30% of the value as- 
sumed do not fit the observed spectra. For small changes 
in the radii the best solutions are always in the 12-14K 
temperature range. 

The temperature derived here is only slightly higher 
than the temperature derived from the SO and CH3OH 
population diagrams. This implies that the observed SO 
and CH3OH transitions are close to thermalization for the 
derived volume density, 2.8 x 10 5 cm -3 . 

3.2. Excitation Temperature and Column Densities 

From the previous analyses the derived temperature of 
the cloud is ~ 13 K. Yet, the rotational transitions of 
the observed molecules may have excitation temperatures 
significantly lower than this value because of subthcrmal- 
ization, i. e., the critical density of the transition is higher 
than the volume density of the cloud, 2.8 x 10 5 cm" 3 . 
In particular this is the case for the high dipole moment 
molecules like HCO+ (3.90 Debyes), HCN (2.98 Debyes), 
H 2 CO (2.33 Debyes) or CS (1.98 Debyes). Thus, for ex- 
ample, the critical density of the H2CO 3o,3-2o.2 and CS 
5-4 are 6 10 6 and 7 10 6 cm" 3 , respectively (Choi & Zhou 



1997) 



If at least two transitions are detected and the emis- 
sion is optically thin, the excitation temperature can be 
obtained from the following two expressions: 



Tj_ _ fi J Vi (T ex ) - J Vi (Tjpg) 1 - e r ' 
Tj fj J Vj (T ex ) - J Vj (T bg ) 1 - e- T i 



and 



9Ki9ii e 



Tj guj gij e hu j/ kT °> 



-(Euj-Eui)/kT e , 



1 



(1) 



(2) 



where the subindices i and j refer to the rotational transi- 
tions, / is the filling factor, J V (T) is the Planck function in 
temperature units, <7k and gj are the K-level and nuclear 
spin degeneracies, respectively, Eu is the energy of the 
upper level of the transition, and r the optical depth. In 
most of the cases here, the spectra of the transitions used 
to derive T cx were obtained at the same angular resolution, 
so fi/ fj — 1. HCO + , and CO are the only molecules with 
rarer isotopes observed. Therefore, for the other molecules, 
their emission is assumed to be optically thin. 



The beam-averaged column densities of the observed 
species were derived assuming that the excitation temper- 
ature, T ox , is the same for all the rotational transitions 
of the same species, with the exception of the SO and 
CH3OH, whose column densities were derived from the 



Population Diagram (§ 3.1.1) 



3.2.1. HCO+ 

The HCO + column density can be derived from the model 
ca rried o ut using the Monte Carlo radiative transfer code 
(§ 3.1.2 ) or directly from the line intensity ratios as de- 
scribed earlier (Eq. || and |l|) . From the model derived in 
§ 3.1.2, the HCO + beam-averaged column density for the 



best fitted model is 1.010 cm . On the other hand, 
for the intensity ratio method, the filling factor ratio of 
the 4-3 and 3-2 transitions, fw/fiW'i should be de- 
rived in order to properly use Eq. ^ and [l| From the 
FWHM of the HCO+ 1-0 emission, 34" (see § |T|) we 
can obtain a rough estimation of this ratio: /21" //2s" — 
(21 2 + 34 2 )/(28 2 + 34 2 ) = 0.8. Thus from the line ratio 
analysis we estimate that the H 13 CO + 3-2 transition is op- 
tically thin (r < 1) and that the excitation temperature is 
within the 6.2-7.0 K range. Therefore, from the H 13 CO + 
3-2 integrated line intensity, the HCO + beam-averaged 
column density is within the 3-5xl0 13 cm" 2 range. The 
line intensity ratio analysis derives an excitation tempera- 
ture lower than the kinetic temperature obtained in § |3.1.2 
and also a lower column density (by a factor 2-3). Uncer- 
taintic s in th e fey //2s" derived from the HCO + 1-0 maps 
(see § 3.1.2| ) could account for this discrepancy. Alterna- 
tively, non-LTE effects, due to subthcrmalization, could 
also cause this discrepancy. 

Since the DCO + has a similar dipole moment as the 
HCO + , we adopted the same excitation temperature de- 
rived from the line ratios of the HCO + : within the 6-7 K 
range, the column density of the DCO + column density 
is 1.0 x 10 12 cm -2 . This value is two orders of magnitude 
lower than that of the HCO + . 



3.2.2. HCN 

For the HCN we fitted the three hyperfine components of 
the 1-0 transitions using CLASS. The best solutions of 
the fit give moderate or low optical depths, r ^ 1.0, al- 
though the excitation temperature is not well constrained, 
T ex > 4 K. However, the observed relative intensities of 
the F=l-1, 2-1 and 0-1 transitions, 1.7±0.2:3.9±0.4:1.0, 
do not agree with the 3:5:1 value expected under LTE 
conditions, but it is within the range of values mea- 
sured in other molecular clouds (Guilloteau & Baudry 
1981| ; Wannier et al. |1974[ Walmsley et al. |1982| ). Veloc- 
ity gradients and a non-uniform density could cause this 
"anomaly" (Guilloteau & Baud ry |1981| ; Cao et al. [1993 ; 
Gonzalez-Alfonso & Cernicharo 1993 ). The 3-2 to 1-0 in- 



tegrated intensity ratio yields an excitation temperature 
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Table 3. Column Densities and Abundances 



Electric Dipole Moment (Debyes) 

Fig. 5. Excitation temperatures as a function of the 
dipole moment for the molecules with two or more transi- 
tions detected towards the quiescent ambient cloud ahead 
of HH 2. The T cx value used for CO, 13 K, is based on the 
assumption that it is thermalized and equal to the kinetic 
temperature derived for the HCO + . 



of ~ 6 K, significantly lower than that derived from the 
population diagram of the SO and CH3OH. This implies 
that the line emission is subthermalized, which could also 
cause the "anomaly" in the 1-0 hyperfine line ratios. The 
HCN beam-averaged column density was derived using 
T cx = 6 K. 

For the DCN, we adopted the same excitation temper- 
ature as the one for the HCN, 6 K. 

3.2.3. H 2 CO 

The two transitions detected yield an excitation temper- 
ature of ~ 6 K, implying that the emission is subthermal- 
ized (the upper limits of the undetected transitions are not 
tight enough to better constrain this value) . As in the case 
of HCN and HCO + this is due to its high dipole moment. 
In order to calculate the H 2 CO beam-averaged column 
density we used this excitation temperature and adopted 
an ortho to para ratio of 1.5, which should be reasonable 
for the H2CO formation and/or ortho-para conversion on 
grain surfaces (Dickens & Irvine 1999). 



3.2.4. CS 

The excitation temperature derived from the two transi- 
tions is ~ 8 K. This slightly higher excitation temperature 
than that derived from the other molecules is possibly due 
to its slightly lower dipole moment. As in the previous 
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cases, this excitation temperature was used to derive the 
column density. 

3.2.5. S0 2 

Of the five SO2 transitions observed, only two transitions 
were detected, 3i,3-2 ,2 and 6o,6-5i,5. If their emission is 
optically thin, then their line ratio gives an unusually high 
excitation temperature, ~ 280 K, much higher than that 
derived from other molecules. At this temperature the 
101,9-92,8 and 10i,g-10o,io transitions should have been 
detected with BIMA. Therefore, this temperature seems 
unreliable. If an excitation temperature of about 10 K 
is assumed, then the line ratio gives an 31,3-20,2 optical 
depth of 10. From the radiative transfer equation for the 
optically thick case, T mD = / (J[T CX ] — J[Tb g ]), the filling- 
factor can be derived and, therefore, a rough estimation 
of the source size, 0, f = 6 2 /(6 2 + 30 2 ) (30" is the angu- 
lar resolution of the spectra). The peak intensity of the 
3i, 3-2o, 2 line, ~ 0.2 K, would imply a source size of only 
5" if the line is optically thick. However, the BIMA maps 
of this transition show that the emission arises from a sig- 
nificantly larger area, ~ 19" (Paper II). Thus, we suggest 
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that the 6o,6~5i,5 line emission is "anomalous", although 
we cannot discern if this is an instrumental effect or if it 
is a real anomaly. In order to estimate the SO2 column 
density we used the 3^3-20.2 transition, assuming an ex- 
citation temperature of 10 K. 

3.2.6. Other molecules 

The CO and its isotopes have a low dipole moment, 0.11 
Debyes, so their emission is likely thermalized. Therefore, 
the excitation temperature used to derived the CO col- 
umn density is 13 K. Figure || shows the values of the 
excitation temperature versus the dipole moment of the 
analyzed aforementioned molecules. With the exception 
of the HCO + , there is a trend of decreasing T cx with in- 
creasing dipole moment, /i, which is expected since higher 
[i implies more severe subthermalization. Fig. |E| can be 
used to roughly extrapolate the excitation temperature 
for the other observed molecules. Thus, we adopted a 
T C x = 10 K for those molecules with /i < 2 Debyes (CN, 
OCS, H 2 S, HCS+, H 2 CS, HNCO, HCOOH and CH 3 SH) 
and T cx = 6 K for the rest of the molecules (0C3H2, 
HC 3 N, CH 3 CN and NH 2 CN), which indeed have a dipole 
moment higher than 3 Debyes. 

Table || shows the beam-averaged column densities de- 
rived from the observations. The fractional abundances of 
the observed molecules with respect to the CO molecule 
are also listed in this table (note that the estimated frac- 
tional abundance relative to the CO is a beam-averaged 
value) . 



4. Discussion 

4-1- Chemical Comparison with other dense molecular 
regions 

In order to be able to establish the effects of the shock- 
induced radiation from HH 2 on the chemical composi- 
tion of the quiescent clump ahead of the HH object, it 
is necessary to compare the relative abundances of the 
molecules detected with respect to other environments as- 
sociated with dense molecular material. Before the radi- 
ation reached the clump, this molecular cloud had likely 
the properties of the starless dark molecular cores. It is 
well known that there are important chemical differences 
between starless dark clouds and star forming regions, and 
that there are important gradients within the same molec- 
ular cloud, a consequence of the evolutionary effects or 
due to the presence of nearby star formation sites (e.g., 
van Dishoeck & Blake 1998). Since we want to compare 



the general chemical properties of the molecular clouds 
with those from the molecular emission ahead of HH 2, 
we have estimated the logarithmic median and standard 
deviation of the relative abundance (with respect to CO) 
of the molecules detected in HH 2 and compare this with a 
sample of three different types of molecular environments: 
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Fig. 6. Relative molecular abundances (with respect to 
the CO) for HH 2 (filled circles), high mass star forming 
cores, low mass star forming cores and dark clouds (i.e., 
without star formation activity). The filled square within 
the bar and the bar show the logarithmic median and the 
logarithmic standard deviation, respectively, of the sam- 
ple. The sample for the high mass star forming clouds is: 
W3 (Helmich fc van Dishoeck |1997|) , M17 and Cepheus 
A (Bergin et al. 1997 ). The sample for the low mass star 



for ming molecular clouds is: IRAS 16293 (v an Di shoeck et 
al. |1995D , NGC 1333 IRAS 4A (B lake et al. |l995D ,the Ser- 
pens S68 (McM ullin et al. |2000| ), 05338-0624 (McMullin, 
Mundy & Blake 1994 ) and several other regions for H2CO 



(Dickens & Irvine 1999) and CH 3 OH (Kalenskii & Sobolev 



1994). The sample for the quiescent dark molecular clouds 



is: OMC-1N (Ungerechts et al. |1997[) TMC-1 (Pratap et 
al. |1997[), L134N (Dickens et al. BOOT), and the CB clouds 



for H 2 CO (Turner 1994). The references of several of the 
selected molecular clouds give the relative chemical abun- 
dances in different positions of the cloud, so most of the 
data plotted in this figure were obtained averaging four 
or more points (the S0 2 and CN values in low mass star 
forming cores are from only one region). 



high mass and low mass star forming molecular clouds, 
and starless dark molecular clouds. Figure ^ shows shows 
this comparison. 

Interestingly, from Figure ^, the observed molecules 
ahead of HH 2 can be classified into four groups, de- 
pending on their abundance with respect to other molec- 
ular clouds. (1) Strongly enhanced molecules: HCO + , 
CH3OH and H 2 CO. Their relative abundances are en- 
hanced mostly within a factor of 5 to 10 with respect to 
the typical values found in the molecular clouds. In partic- 
ular, the enhancement of HCO + and H2CO in HH 2 with 
respect to the high mass star forming regions is more than 
an order of magnitude, but only by a factor of 2-3 with re- 
spect to the dark molecular clouds. (2) Weakly enhanced 
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molecules, SO and SO2. Their enhancement with respect 
to the molecular clouds is a factor of 2-4. The only ex- 
ception is that the SO2 has similar relative abundances as 
the massive star forming regions. (3) The CN shows no 
apparent enhancement with respect to the dark molecu- 
lar clouds, whereas it is depleted with respect to the high 
mass star forming regions. The strong enhancement (an 
order of magnitude) with respect to the low mass star 
forming regions may not be significant since there is only 
one data point. (4) Depleted molecules: CS and HCN. In 
HH 2 these molecules are clearly depleted with respect to 
the other molecular clouds, with depletion factors in the 
3 to 7 range. 

4-2. Comparison with the VW99 Models 

Here we briefly compare our observational results with the 
abundances derived in the VW99 models. A more detailed 
model of the HH2 clump (in light of the present observa- 
tions) will be the subject of a future paper. Columns 3 
to 6 of Table ^ lists the best matched theoretical column 
densities and the visual extinction at which this column 
density was obtained for two models out of the grid from 
VW99 for clump densities of 1 10 4 cm" 3 (Model 4) and 
of 1 10 5 cm -3 (Model 5), 3 years after irradiation of the 
clump has started. 

Table [| clearly confirms the qualitative results from 
VW99 where many species were predicted to be abundant 
in clumps ahead of HH objects. The VW99 models were 
not specific to any particular HH object but use an exten- 
sive gas and grain phase chemistry to follow the formation 
of a clump ahead of a shock with a fixed enhanced radia- 
tion field (equivalent to 20 times the ambient interstellar 
radiation field) acting upon it. The clump is treated as a 
one-dimensional slab extending up to 6 visual magnitudes 
to the center. 

In general, the VW99 models suggest that: (i) HH2 is 
a young object, probably on the order of 1000 years; (ii) 
the best estimate for the volume density of the clump is 
< 10 5 cm -3 , somewhat less than the density derived from 
observations; (iii) the visual extinction at the center of the 
SO peak emission with respect to HH 2 is about 5-7 mags; 



VW99 models only extend to 6 mag; although quantita- 
tively the models do not match the observed abundances 
for every species at the maximum Ay, most of them are 
well represented at some visual extinction. In fact, Model 4 
of VW99 reproduces the observed abundance (within half 
order of magnitude) for every species for an optimal visual 
extinction, different from species to species: this may be an 
indication of the displacement found with the BIMA array 
(Paper II). More specifically, it is interesting to note that 
HCO + , the best tracer of HH clumps, is underabundant 
for visual extinction larger than ~ 1 mag. This is because 
at high visual extinction the clump is more shielded from 
the radiation and this would limit the amount of ionized 
carbon, yielding a low production of HCO + ; a higher radi- 
ation field may well therefore increase its abundanc e. The 
presence of X-ray emission in HH 2 (Pravdo et al. 2001) 
may have some effects on the chemistry. A detailed model 
for HH 2 which takes such considerations into account is 
in preparation. 



5. Conclusions 

We presented a 3-0.8 mm molecular line survey toward the 
quiescent molecular clump ahead of the bright HH 2 ob- 
ject. A total of 14 species were detected (including differ- 
ent isotopes and deuterated species) and the upper limits 
of 12 more species were obtained. Multi-transitions ob- 
servations of HCO + , SO and CH3OH show that the gas 
is not only quiescent, but also cold, ~ 13 K, in spite of 
being close to HH 2. The density of the cloud is roughly 
3 10 5 cm" 3 . Comparisons of the relative abundances (with 
respect to CO) were made with different dense molecu- 
lar gas environments, including quiescent, starless dark 
clouds, and active low and high mass star forming cores. 
We confirm the peculiar chemical composition of the gas 
ahead of the HH objects, a consequence of the mantle re- 
moval and further photochemistry produced by the strong 
radiation generated in the shocks. In particular, we found 
that: 

— HCO + , CH3OH and H2CO are strongly enhanced with 
respect to other molecular cores. 

— SO and S0 2 are weakly enhanced with respect to other 
molecular environments. The only exception is that the 
SO2 has similar relative abundances as the massive 
star forming regions. 

— CN does not show a significant enhancement or deple- 
tion with respect to the starless dark clouds. However 
it is underabundant with respect to high mass star 
forming regions. 

— HCN and CS show clear depletion when the HH 2 rela- 
tive abundances are compared with most of the molec- 
ular clouds. 

Finally, the chemical composition of HH 2 confirms the 
qualitative results of the VW99 complex chemical model 
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that follows the chemical behavior of a molecular clump 
irradiated by a HH object. 
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Appendix A: The Population Diagram Analysis 

The population diagram (also called rotation diagram in 
the literature) analysis has been extensively used in the 
literature to obtain the temperature and the column den- 
sities for molecules with several transitions observed (e. g., 
Linke, Frerking & Thaddeus [l979| ; Turner [I991D . In the 
use of the population diagram, we assume LTE condi- 
tions (i.e., all the transitions are populated with a single 
excitation temperature T ox ), optically thin emission and 
T C x 3> Tb g . Under these assumptions the population dia- 
gram can be described with the following expression: 



log / = log 



N oa _ loge 

Qrnt T 



E,, 



with 



I = 



3fc / T mh dv 
8Tr 3 g K giSn 2 is 



(A.2) 



where gK and g\ are the K-level and nuclear spin degen- 
eracies respectively, \i the electric dipole moment of the 
molecule, S, u, T mD and E u are the line strength, fre- 
quency, main-beam brightness temperature and the en- 
ergy of the upper state (in temperature units) of the 
transition, N\, a the total beam-averaged column density 
of the molecule, Qrot the rotational partition function of 
the molecule. This formula appears in the population di- 
agram (in which the left hand part of the formula is plot- 
ted versus E u ) as a straight line. However, several factors 
can yield temperature and column densities that are sig- 
nificantly different from their true values: subthermaliza- 
tion (or non-LTE excitation), significant optical depths, 
excitation gradients within the beam and beam dilution 
(when diff erent angular resolutions are used). Goldsmith 
& Langer (|l999|) describe the effects that subthcrmaliza- 
tion and non-null optical depths have on the population 
diagram. 

In order to better use the population diagram tech- 
nique we just assume that all the rotational levels are 
populated with the same excitation temperature (which 
in most cases is likely lower than the kinetic temperature: 



see Goldsmith & Langer 1999). In this case the expression 
is: 



log / = log ■ 



N h 



loge 
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E u -\ogC r -\ogF T 



(A.3) 
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Fig. A.l. Ft as a function of temperature for three dif- 
ferent frequencies. 



(A.l) C T 



with C T and Ft are defined as: 
r 



l-e~ 



and 



Ft 



J v [T c: 



(A.4) 



(A.5) 



J„[T] is the Planck function in temperature units. Note 



that when r <C 1 and T cx 3> T,„, equation A.3 becomes 



equation A.l. However, for those regions where the tem- 
perature is low and with moderate densities, Ft should be 
taken into account, especially for the low frequency transi- 
tions (see Fig. A.l). C T should also be taken into account 
even for moderately small opacities: e.g., C T — 1.25 for 
r = 0.5. 

The optical depth for each transition can be estimated 
if true column density, Ntrue, is known. Since N tm c X / = 
N oa , the optical depth of a given rotational transition is 
given by: 



t .j,j+i = 



87r 3 /i 2 A^ ba J e - aJ ( J + 1 )(e 2aJ - 1) 

3h f Qrot 



(A.6) 



Therefore, if the filling factor / is included as an additional 
free parameter (with T and N Da ) , the optical depth can be 
corrected whenever using the population diagram method. 

To estimate the best set of solutions we calculated the 
X 2 value for each set of parameters, T YOt , N D& and /, taking 
into account the uncertainties (rms noise and calibration 
errors) in a similar way as that described by Nummelin et 
al. ( p998| ) and Gibb et al. ( p00C| ): 



* 2 = E 

i=i 



robs T> 

_z _i 

_obs 



cak 



(A.7) 
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For those transitions undetected, the 3-er upper limit for 
each transition was used whenever I^ alc < 3 x a. The 
1-cr uncertainties for T rot and TVba were obtained from 
the 68% confidence region 

of x 2 (e.g., Gibb et al. |200(i| ), 
which is enclosed within the Xmm + ^X 2 interval, with 
Ax 2 = x/(N — P) where N is the number of data points 
(including those undetected transitions that satisfy the 
aforementioned condition) , P is the number of free param- 
eters (3 when the opacity correction is taken into account 
and 2 when optically thin emission is assumed), and x 
equals to 2.3 and 3.5 for 2 and 3 fre e para meters, respec- 
tively (Lampton, Margon & Bowyer 1976 ). 
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